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1.0  INTRODUCTION 

Computed  tomography  (CT)  has  exhibited  great  potential  for  the  determination  of 
temperature  and/or  specie  concentration  in  nonaxisymmetric  gas  flows  (Refs.  1-4).  With  the 
advent  of  vectoring  nozzles  in  turbine  engines  and  the  use  of  multiple  rocket  engine  nozzle 
combinations,  a  technique  for  determining  the  spatial  profiles  of  static  temperature  and/or 
specie  concentration  has  become  important.  These  properties  are  required  for  code  validation, 
heat  load  calculations,  and  signature  studies.  Point  property  measurement  techniques  such 
as  laser-induced  fluorescence  also  show  promise  for  application  to  gas  diagnostics  (Refs.  5-8). 
Point  property  techniques  are  advantageous  since  they  do  not  require  the  determination  of 
point  properties  indirectly  from  line-of-sight  measurements.  However,  the  complexity  of  the 
equipment  needed  for  point  property  techniques  such  as  laser-induced  fluorescence  is  far 
greater  than  that  required  for  emission  tomography.  Likewise,  absorption  tomography  also 
requires  a  substantial  investment  in  equipment.  Both  these  techniques  require  expensive  and 
complicated  tunable  laser  sources  and  their  associated  beam  delivery  systems.  Computed 
emission  tomography,  on  the  other  hand,  can  be  applied  with  measurement  technologies 
presently  available  such  as  low-light-level  cameras  or  spatially  scanning  bandpass-filtered 
radiometers. 

Presented  here  is  a  summary  of  an  investigation  into  the  applicability  of  emission 
tomography  for  gas  diagnostics.  The  goal  of  this  study  was  to  determine  if  existing  computed 
tomography  algorithms  could  be  employed  to  reconstruct  emission  cross-sectional  profiles 
from  emission  data  on  high-temperature  gases,  and  whether  these  data  could  be  utilized  to 
yield  temperature  and/or  specie  partial  pressure  profiles. 

2.0  COMPUTED  TOMOGRAPHY 

Computed  tomography,  in  a  general  sense,  is  a  mathematical  technique  for  computing 
a  point  property  from  line-of-sight  integrals  of  the  same  property.  Consider  Fig.  1,  in  which 
the  line  integral  of  the  property  c  is  determined  for  several  different  paths  denoted  by  the 
angles  0  and  a.  The  path  integrals  of  c,  denoted  c(0,  a)  are  given  by 

<■) 

where  P(0,  a)  is  the  path  determined  by  the  angles  6  and  a.  These  values  serve  as  the  input 
to  a  tomographic  reconstruction  algorithm  which  would  generate  point  properties  of  c,  i.e., 
c  =  c(x,y).  Therefore,  tomography  can  be  thought  of  as  the  reverse  process  of  Eq.  (1)  in 
which  c(x,y)  is  found  from  c(0,  o).  t^o  a  priori  knowledge  of  the  function  c(x,y)  is  required. 
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The  most  readily  accessible  set  of  tomographic  reconstruction  algorithms  is  the  Conner 
algorithms  available  from  the  Conner  laboratory  of  Lawrence  Berkeley  Laboratory  (Ref. 
9).  Although  these  algorithms  were  written  primarily  for  the  medical  community,  they  have 
been  employed  in  other  areas  such  as  nondestructive  testing.  The  theory  of  tomographic 
reconstruction  is  not  presented  here,  but  there  are  various  numerical  approaches  to  performing 
Radon  transforms  such  as  direct  algebraic  inversion  techniques,  two-dimensional  Fourier 
transform  analysis,  and  iterative  techniques.  The  Conner  computer  code  contains  all  the  major 
reconstruction  techniques,  as  well  as  means  to  handle  several  different  source  and/or  detector 
geometries. 


3.0  RADIATIVE  TRANSFER 

An  optical  gas  diagnostic  technique  employing  CT  must  satisfy  two  basic  requirements. 
One  requirement  is  that  the  reconstructed  quantity  must  be  a  function  of  temperature,  number 
density,  or  other  desired  flow  parameters.  The  other  requirement,  which  is  essential  for  the 
application  of  tomography,  is  that  from  the  measured  quantity,  a  line-of-sight  integral  of 
the  reconstructed  quantity  must  be  obtained  which  will  serve  as  input  for  a  tomographic 
reconstruction.  For  example,  consider  a  measurement  of  transmittance.  The  transmittance 
is  not  a  line-of-sight  integral  of  some  quantity.  However,  a  line-of-sight  (LOS)  integral  of 
the  absorption  coefficient  can  be  extracted  from  the  transmittance  measurements  utilizing 
Beer’s  law; 


7(0, «)  =  e-K(e.“>  =  (2) 

where  r  is  the  transmittance,  ka(x,y)  is  the  absorption  coefrident,  K(6,  a)  is  the  LOS- 
integrated  absorption  coefficient,  and  P(6,a)  represents  the  integration  path  as  shown  in 
Fig.  1.  The  integrated  absorption  coefrident  is  easily  calculated  from  the  measured 
transmittance: 


=  -hi7(e,a)  (3) 

Since  this  is  an  LOS  integral,  it  is  applicable  to  tomographic  reconstruction,  and  since  the 
absorption  coeffident  is  temperature-  and  number  density-dependent,  the  transmittance  is 
a  possible  tomography  gas  diagnostic  tool  (Refs.  10-11).  To  investigate  the  applicability  of 
emission  measurements  for  tomography  reconstruction,  consider  radiative  transfer  through 
a  nonhomogeneous  gas.  The  radiative  transfer  equation  is  given  by: 

N(X)  =  +  No(X)je "  (4) 
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where  the  symbols  are  defined  as  follows: 


Aif  Transition  probability,  sec-' 

c  Speed  of  light,  cm-sec“' 

h  Planck's  constant,  erg-sec 

Oj  Number  density  of  emitting  atoms  in  initial  state  i,  cm-^ 

X  Wavelength  of  spectral  line,  nm 

N(X)  Emitted  spectral  radiance,  w/(cm^-sr-nm) 

No(X)  Incident  spectral  radiance,  w/(cm2-sr-nm) 

ka  Absorption  coefficient,  cm-' 

f(X)  Spectral  line  shape  function,  nm-' 

1  LOS  path  length 


The  absorption  coefflcient  above  includes  the  stimulated  emission  term  and  is  given  by: 


nA4lO-2if(X)Air  ,  ,  . 


(5) 


where  nf  is  the  number  density  of  atoms  in  the  final  (lower)  state  and  T  is  the  temperature. 
Assuming  the  gas  is  in  thermal  equilibrium,  the  number  density  of  emitting  atoms  in  state 
i  is  distributed  according  to  a  Maxwell-Boltzmann  distribution  and  is  given  by: 

g.g  -  Ejhc/kT 


where  the  symbols  are  defined  by: 

Ej  Energy  of  ith  state,  cm-' 

Q(T)  Partition  function 

T  Temperature,  K 

k  Boltzmann’s  constant,  erg  K-' 

n  Total  number  density  of  emitting  atoms,  cm“^ 

gi  Degeneracy  of  i'**  state 


The  spatial  dependence  of  the  spectral  radiance  is  through  the  temperature  and  number  density 
which  are  position  dependent.  Examination  of  Eq.  (4)  indicates  that  from  a  measurement 
of  the  emitted  spectral  radiance  alone,  one  cannot  determine  a  simple  line  integrable  quantity. 
However,  if  it  is  assumed  the  gas  is  optically  thin  (i.e.,  no  self  absorption,  k,  =  0)  and  there 
is  no  incident  radiation  (N^  =  0),  Eq.  (4)  reduces  to: 


N(X)  = 


AifhclO’f(X)  rf 


(7) 
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If  one  assumes  that  the  spectral  instrument  response  function,  g(X),  of  the  detecting  system 
is  wide  in  comparison  to  the  spectral  line  width  of  the  emitting  radiation,  the  measured  radiance 
is  approximately  given  by: 


AjfhclO’g(Xo)  ft 

N  =  - — - 1  nrdx 

4ffXo  Jo  ‘ 


(8) 


where  Xo  is  the  emitting  spectral  line  position. 

A  comparison  of  Eq.  (8)  and  Eq.  (1)  indicates  the  measured  radiance  has  the  desired  LOS 
integrability  for  use  as  input  to  computed  tomography.  Note  the  reconstructed  quantity  is 
a  multiplicative  factor  of  the  number  density  of  the  radiating  atoms  in  the  initial  quantum 
state.  Equation  (6)  expresses  this  as  the  total  number  density  of  the  emitting  atoms  and  the 
temperature.  Therefore,  it  is  possible  to  determine  the  number  density  of  the  radiating  species 
from  the  reconstructed  state  specific  number  densities,  provided  the  spatial  profile  of  the 
temperature  is  known. 


Measurement  of  two  radiances  originating  from  two  different  spectral  lines  of  the  same 
species  could  be  utilized  to  simultaneously  determine  both  number  density  and  temperature 
under  certain  instances.  Consider  two  spectral  lines,  X^  and  X^,  from  the  same  specie 
originating  from  two  different  energy  levels,  E;  and  Ei'-  Reconstruction  of  the  measured 
radiances  would  yield  the  following  reconstructed  quantities: 


c(x,y) 


AifhclO’g(Xo) 

4tXo 


ni(x,y) 


(9) 


and 


c'(x,y) 


Aj'fhclO^giXb) 

4tX6 


ni(x,y) 


(10) 


Taking  the  ratio  of  the  two  reconstructed  quantities  yields: 

c(x,y)  ^  Aifg(Xo)X6  nj(x,y) 
c'(x,y)  Aj'fg(X6)Xo  n'i(x,y) 

which,  with  the  help  of  Eq.  (6),  can  be  rewritten  as 

e-(Ei  -  EiVkT 

c'(x,y)  Aifg(X6)Xog'i 


(11) 


(12) 


Values  for  the  Einstein  spontaneous  emission  coefficient,  A,  degeneracy,  and  wavelengths 
can  be  found  in  the  literature,  and  the  instrument  function  can  be  established  from  the 
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instrument  calibration  (Ref.  12).  Once  these  parameters  are  established,  the  temperature 
equation  [Eq.  (12)]  above  can  be  solved  for  temperature: 


T(x,y) 


[c(x.y)XoAif8(X6)gi] 

lc'(x,y)XoAjfg(Xo)gjll 


(13) 


Once  the  temperature  is  known  at  each  spatial  location,  Eqs.  (6)  and  (9)  (or  10)  can  be  employed 
to  determine  the  number  density:  s 


n(x,y) 


c(x.y)4>rX^^Q(T) 
Ajfhc  10^g(Xo)gje  ~ 


(14) 


4.0  EXPERIMENTAL  APPROACH 


The  sodium  doublet  at  S89  and  589.  S  nm  was  chosen  for  the  emission  tomography 
measurements.  This  doublet  appears  as  a  single  spectral  line  at  moderate  resolution  and  has 
been  observed  to  occur  in  varying  degrees  in  most  engine  plumes.  The  radiation  from  sodium 
at  this  wavelength  and  from  the  potassium  doublet  at  766.5  and  769.9  nm  is  consistently 
observed  in  the  hydrogen-oxygen  space  shuttle  main  engine  (SSME)  plume.  This  flame  cannot 
be  duplicated  in  the  lab,  but  a  hydrogen-air  flame  can  be  created  utilizing  a  small-scale  burner. 
Since  both  sodium  and  potassium  are  excellent  radiators,  the  number  density  required  to 
yield  the  observed  radiation  levels  in  the  SSME  is  extremely  low.  An  analysis  of  SSME  signature 
data  indicated  that  the  radiance  due  to  sodium  and  potassium  was  approximately  equal  at 
the  nozzle  exit  plane;  however,  the  concentration  required  to  yield  the  observed  radiation 
level  was  determined  to  be  0.87  and  0.06  ppb,  respeaively  (Refs  13).  Recent  measurements 
on  MMH-N2O4  (monomethylhydrazine)  and  UDMH-N2O4  (unsymmetrical  dimethylhydra- 
zine)  propellant  engines  have  indicated  that  they  also  possess  sodium  emissions  in  their  plumes. 

The  589-  to  S89.5-nm  doublet  is  the  only  spectral  feature  of  sodium  or  potassium  that 
is  consistently  and  easily  observed  at  moderate  spectral  resolution  in  exhaust  plumes.  This 
is  unfortunate,  since  it  precludes  the  use  of  emission  measurements  from  two  spectrally  resolved 
lines  of  the  same  gas  to  simultaneously  determine  both  temperature  [from  Eq.  (13)]  and  number 
density  [from  Eq.  (14)].  Therefore,  a  temperature  profile  must  be  determined  independently. 
This  temperature  profile  can  then  be  employed  to  determine  the  number  density  from  the 
reconstructed  emission  coefficients  by  using  Eq.  (14). 

A  proof-of-principte  experiment  was  designed  around  the  use  of  a  flat  flame  burner.  The 
particular  burner  used  was  a  2-in.-square  burner  shown  in  Fig.  2  obtained  from  Research 
Technologies.  The  body  of  the  burner  consisted  of  a  Hastalloy®  honeycomb  interlaced  with 
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0.02^in,-lD  hypodermic  tubes  shown  in  Fig.  3.  The  hydrogen  was  delivered  by  the  hypodermic 
tubes,  and  the  oxidizer  by  the  voids  in  the  honeycomb.  Hydrogen  and  air  flow  rates  were 
measured  using  Edwards-Datametrics  50  and  100  slpm  (standard  liters  per  minute)  full-scale 
digital  flow  meters,  respectively. 

4.1  SODIUM  LASER  ABSORPTION  MEASUREMENTS 

The  first  measurements  on  the  burner  system  were  sodium  spectral  absorption 
measurements  to  determine  the  extent  of  self  absorption.  The  configuration  for  this  experiment 
is  shown  in  Fig.  4.  The  source  utilized  is  an  argon-ion  laser-pumped  tunable  ring  dye  laser. 
The  dye  laser  output  was  split  between  a  laser  fluorescence  reference  cell  and  the  probe  beam 
which  traversed  the  fiat  flame  burner  approximately  0.75  in.  above  the  burner  base.  The 
signal  from  the  fluorescence  reference  cell  was  used  to  indicate  when  the  laser  had  been  tuned 
to  the  sodium  absorption  line.  These  measurements  were  conducted  over  a  wide  range  of 
flow  conditions  and  no  measurable  absorption  was  found.  This  implies  that  the  necessary 
criterion  for  the  use  of  emission  tomography  is  met;  that  is,  an  optically  thin  medium  is 
assumed. 

4.2  TOMOGRAPHY  EMISSION  MEASUREMENTS 

Tomographic  emission  data  were  acquired  using  a  low-light-level  camera.  The  camera 
selected  was  a  Burle  model  07985  since  it  is  equipped  with  an  automatic  gain  control  (AGC) 
enable/disabie  option.  This  feature  allows  acquisition  of  images  under  varying  light  conditions 
that  are  indicative  of  the  true  light  level.  The  camera  was  equipped  with  a  Nikon  f/4  80-2(X) 
mm  zoom  lens  which  gave  an  image  of  the  burner  which  filled  approximately  75  percent 
of  the  video  frame  at  a  source  distance  of  51  in.  Mounted  to  the  front  of  the  lens  was  a 
bandpass  optical  filter.  The  transmission  curve  of  the  filter  is  shown  in  Fig.  5  from  which 
the  center  wavelength  and  bandpass  are  found  to  be  588  and  5  nm,  respectively.  The  output 
of  the  low-light-level  camera  was  recorded  on  VHS  video  tape  which  was  later  played  back 
for  analysis. 

To  obtain  images  of  the  flame  at  various  angles,  the  flat  flame  was  mounted  on  a  rotation- 
translation  stage  as  shown  in  Fig.  2.  This  allowed  the  burner  to  be  rotated  through  180  deg. 
Figure  6  is  a  schematic  diagram  of  the  experimental  setup,  and  the  photograph  in  Fig.  2  shows 
the  position  of  the  camera  in  the  upper  left  and  its  relation  to  the  burner.  All  data  were  taken 
0.75  in.  above  the  burner  face. 

Processing  of  the  video  was  performed  on  a  Quantex  image  processor.  Individual  frames 
corresponding  to  the  various  angles  were  digitized.  From  these  frames,  an  individual  horizontal 
scan  corresponding  to  a  slice  0.75  in.  above  the  burner  face  was  taken,  as  illustrated  in  Fig. 
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6.  These  data  then  served  as  input  to  the  tomography  algorithms.  Data  were  taken  every 
12  deg  over  180  deg,  with  each  slice  consisting  of  233  pixels.  The  flow  rates  through  the  burner 
were  16.37  slpm  hydrogen  and  36.47  slpm  air,  which  yields  an  equivalence  ratio  of  l.Jl. 
Flow  settings  were  chosen  to  achieve  the  maximum  sodium  emission;  this  was  by 

observing  the  camera  video  output.  Repeated  operation  of  the  burner  indicated  that  maximum 
radiation  in  the  S88-nm  region  was  not  consistent  with  the  flow  settings. 

The  horizontal  slices  of  the  flame  image  were  used  as  input  to  the  Donner  tomographic 
reconstruction  algorithms.  The  fan-beam  geometry  option  of  the  Donner  algorithms  was 
selected  for  data  reduction  and  the  conjugate  gradient  technique  for  reconstruction.  Appendix 
A  is  a  listing  of  the  main  subroutine  used  to  invoke  the  tomography  algorithms.  Figure  8 
illustrates  the  resulting  reconstructed  itnages  which  can  be  compared  to  the  color  photograph 
(Fig.  7)  of  the  Hj-air  flame.  Note  the  nonuniformity  and  strong  emission  characteristics  in 
both  the  photograph  and  reconstructed  image  which  appears  in  the  periphery  of  the  flame. 

Contrary  to  our  initial  assumptions,  the  flame  was  far  from  uniform.  This  discovery  poses 
a  substantial  problem  from  the  standpoint  of  modeling.  If  the  fuel-air  mixture  is  not  uniformly 
distributed  across  the  burner  face,  then  using  the  adiabatic  flame  temperature  from  the  NASA 
273  equilibrium  chemistry  code  for  the  flame  temperature  would  be  in  error  (Ref.  14).  An 
attempt  was  made  to  measure  the  flame  temperature  by  traversing  a  Pt-Pt/Rh  thermocouple 
through  the  flame.  Figure  9  illustrates  the  results  of  the  thermocouple  traverses.  The 
temperature  shown  in  Fig.  9  is  the  thermocouple  temperature,  not  the  gas  temperature.  A 
difference  exists  between  these  measurements  because  of  the  convective-radiative  heat-transfer 
losses.  Two  features  are  apparent  from  examination  of  the  thermocouple  data:  the  nonuniform 
distribution  and  the  low  measured  temperature  values.  The  measured  temperatures  are  lower 
than  the  adiabatic  flame  temperature,  2,404  K,  by  approximately  1,200  K.  This  difference 
is  far  greater  than  could  be  explained  by  convective-radiative  heat-transfer  losses  from  the 
thermocouple  which  would  introduce  an  error  of  approximately  100  K.  An  explanation  can 
be  found  in  examining  Fig.  7,  Note  the  radiation  from  the  burner  face  (an  extremely  bright 
glow)  which  implies  a  substantial  heat  transfer  to  the  burner  face  and,  therefore,  may  explain 
the  lower  measured  flame  temperature.  Nonuniformities  in  the  flow  were  initially  thought 
to  be  due  to  nonuniform  distributions  in  the  mixing  of  fuel  and  oxidizer. 

To  examine  the  nonuniformity  in  the  fuel-air  ratio,  a  series  of  mass  spectrometer 
measurements  was  made  across  the  burner  face.  A  I/8-in.-diam  probe  attached  to  a  mass 
spectrometer  was  traversed  across  the  burner  0.75  in.  above  the  burner  face.  The  hydrogen 
fuel  was  replaced  by  argon  and  the  air  oxidizer  was  replaced  by  N2.  Figure  10  illustrates  the 
results  of  the  measured  argon/N2  ratio  profiles.  Note  that  in  the  main  center  portion  of  the 
burner  area,  a  uniform  mix  was  obtained.  Data  in  the  periphery  of  the  flow  are  less  accurate 
due  to  the  small  concentration  of  Ar,  where  as  much  as  a  factor  of  two  difference  is  seen 
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in  the  different  regions  of  the  burner.  Based  on  the  mass  spectrometer  profiles,  particularly 
in  the  0-  to  20-mm  region  of  the  burner,  the  temperature  nonunifoimities  cannot  be  attributable 
to  poor  mixing. 

While  investigating  the  sodium  radiation,  researchers  discovered  that  the  burner  face 
temperature  played  an  important  role  in  the  sodium  emission  level.  Further  investigation 
led  to  the  discovery  that  placing  an  object  such  as  a  stainless  steel  tube  in  the  flame  can  induce 
stronger  sodium  emission.  If  the  tube  is  cooled,  the  sodium  radiation  essentially  disappears; 
however,  when  allowed  to  reheat  in  the  flame,  a  bright  sodium  radiation  is  observed.  Figure 
1 1  is  a  spectrum  obtained  by  employing  an  optical  multichannel  analyzer  (OMA)  to  view 
the  region  above  the  heated  tube.  The  spectrum  illustrates  the  presence  of  both  sodium  and 
potassium  in  the  flame  above  the  heated  tube.  When  cooled  by  flowing  water  through  the 
tube,  the  radiation  from  the  sodium  and  potassium  decreased  in  intensity  until  it  was 
indistinguishable  from  the  background  noise,  although,  to  a  lesser  extent,  the  same  result 
was  obtained  from  the  burner  itself.  When  the  burner  was  operated  at  the  highest  flow  rates 
achievable  with  the  present  system,  the  flame  partially  detached  from  the  burner  face.This 
decreased  the  heat  transfer  to  the  burner,  and  a  marked  decrease  in  the  sodium  radiation 
levels  was  measured.  Conversely,  lowering  the  overall  gas  flow  rates  raised  the  burner  face 
temperature,  and  a  substantial  increase  in  sodium  and  potassium  emission  levels  was  observed. 
Although  there  is  no  explanation  of  this  phenomenon  at  present,  the  mechanism  responsible 
for  sodium  excitation  is  important.  If  sodium  is  to  be  used  as  a  diagnostic  tool  in  tomography 
or  laser-induced  fluorescence,  the  understanding  of  its  excitation  mechanism  is  important 
for  use  in  gas  diagnostic  modeling  computer  codes. 

5.0  SUMMARY 

The  radiative  transfer  equations  for  atomic  spectral  line  radiation  emanating  from  a  hot 
gas  were  examined  for  applicability  to  tomography.  Methodologies  for  determining 
temperature  and/or  number  densities  from  the  reconstructed  quantities  were  derived.  A  proof- 
of-principle  experiment  was  performed  to  demonstrate  the  applicability  of  emission 
tomography.  Tomography  data  consisting  of  atomic  sodium  radiation  from  a  H2-air  flame 
from  a  flat  flame  burner  were  acquired  using  an  optically  filtered  low-light-level  camera. 
The  sodium  radiation  data  were  reconstructed  by  application  of  the  Conner  tomography 
algorithms.  Based  on  visible  photographs  and  visual  examination  of  the  flame,  the 
reconstructed  image  was  found  to  be  consistent  with  observations.  However,  large 
nonuniformities  in  the  flame  were  observed.  These  nonuniformities  and  discretions  between 
equilibrium  chemistry  code  results  and  measured  thermocouple  data  precluded  a  quantification 
of  the  sodium  number  density  or  temperature  from  the  data.  Farther  investigation  into  the 
source  of  the  nonuniformities  suggested  that  a  surface  reaction  mechanism  is  responsible 
for  the  excitation  of  sodium  in  the  flame.  At  this  time  the  exact  nature  of  this  mechanism 
is  unclear. 
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6.0  RECOMMENDATION 

A  quantiHable  flame  is  required  to  thoroughly  test  the  applicability  of  tomography  for 
diagnostics.  Although  it  has  been  shown  that  it  is  possible  to  reconstruct  emission  images 
consistent  with  visual  observation,  a  more  detailed  and  controlled  experiment  is  required 
to  test  the  ability  of  emission  tomography  for  temperature  and  number  density  determination. 
The  sodium  excitation  process  appears  to  be  surface  related  and  not  quantifiable  at  this  time. 
Therefore,  it  would  be  advantageous  to  use  another  specie  that  can  be  easily  seeded  into 
the  fuel  or  oxidizer.  If  a  seed  specie  can  be  identified  which  will  equilibrate  to  the  flame 
temperature  and  has  sufficient  emission  characteristics,  a  more  detailed  experimental  study 
would  be  in  order.  This  may  require  the  use  of  a  higher  temperature  flame  such  as  a  hydrogen* 
oxygen  flame  to  create  sufficient  excitation  energies  to  excite  the  seeded  specie(s).  Such  an 
experiment  requires  sufficient  gas  flow  capabilities  to  achieve  gas  velocities  above  the  burner 
face  high  enough  to  lift  the  flame  front  off  the  face  of  the  burner  and  lower  the  heat  transfer 
to  the  burner.  Under  these  circumstances,  a  modeling  of  the  flame  chemistry  and  temperature 
should  be  sufflciently  accurate  to  use  as  a  control  in  an  experiment  to  test  temperature  and/or 
number  density  diagnostics  using  emission  tomography. 
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Figure  1.  Geometry  model  describing  relations  of  tomographic  quantifies. 
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Figure  3.  Photograph  of  flat  flame  burner  face  showing  array  of  hypodermic  tubes  in  Hastalloy® 
honeycomb. 
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Figure  4.  Sodium  laser  absorption  experimental  configuration 
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Figure  5.  Sodium  optical  filter  transmission  curve. 
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Figure  6.  Schematic  of  tomographic  emission  experiment  and  the  digitized  field  of  view 
required  for  tomographic  reconstruction. 
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Kiuiirc  7.  Non-fiUcrctl  plioloiiraph  of  Hi-air  flame  on  flat  flame  hiirner. 
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Figure  8.  Reconstructed  sodium  emission  profile  from  HiO-air  flat  flame  liurner  iitiliT^ing 
bandpass-filtered  camera  data. 
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Figure  10.  Ratio  of  abundance  of  air  to  N2  above  flat  flame  burner. 
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APPENDIX  A 
DONNER  ALGORITHMS 

The  Oonner  algorithms  for  reconstruction  tomography  (Ref.  9)  are  a  library  of  subroutines 
for  three-dimensional  reconstruction  problems  that  arise  in  the  medical  and  physical  sciences. 
These  subroutines  were  written  at  the  Lawrence  Berkeley  Laboratory  of  the  University  of 
California,  Berkeley,  CA.  They  are  written  in  the  FORTRAN  programming  language  and 
are  in  the  public  domain.  The  subroutines  have  been  run  at  AEDC  on  the  AMDAHL 
mainframe  computer  and  the  IBM  RISC  60(X)  workstation.  For  industrial  computed 
tomography  (CT),  the  projection  data  input  into  the  subroutines  are  normally  obtained  from 
transmitted  intensity  profiles  produced  as  X-rays  are  transmitted  through  an  object.  In  addition 
to  transmission  data,  the  medical  field  also  uses  emission  tomography  to  determine  the 
concentrations  of  a  radionuclide  within  the  body.  As  discussed  earlier  in  this  report,  it  is 
not  necessary  for  the  radiation  to  be  X-rays  or  gamma  rays.  The  mathematics  and  algorithms 
apply  to  visible  light  as  well. 

Based  on  earlier  work  here  at  AEDC,  three  of  the  nine  reconstruction  algorithms  in  the 
Donner  library  provided  the  best  reconstructions  for  our  requirements  in  terms  of  speed  and 
memory  requirements.  The  classes  of  algorithms  designated  as  filtered-back  projection  (the 
subroutine  FILBK)  and  convolution  (subroutine  CONVO)  are  most  commonly  used  in 
commercial  CT  systems  because  of  their  high  speeds.  The  image  quality  is  excellent  if  they 
are  given  a  full  complement  of  input  data  with  high  signal-to-noise  ratios.  The  iterative-type 
algorithms  (subroutines  CONOR  and  GRADY)  are  not  used  in  commercial  systems  because 
of  their  long  reconstruction  times.  When  CT  data  are  acquired  under  dynamic  conditions 
such  as  rocket  motor  or  turbine  engine  tests,  It  is  not  feasible  to  acquire  a  full  complement 
of  data,  and  there  is  a  variety  of  possible  sources  of  noise.  The  iterative-type  algorithms  appear 
to  produce  better  images  from  limited  data  sets;  thus,  they  are  used  most  often  at  AEDC 
to  reconstruct  CT  images  from  test  data. 

Listed  on  the  following  pages  is  a  typical  ‘main  program’  used  to  call  a  reconstruction 
routine.  The  main  program  allows  the  user  to  set  the  geometry  parameters  required  by  the 
reconstruction  subroutine. 

Before  any  of  the  reconstruction  routines  are  called,  the  user  must  call  the  subroutine 
SETUP.  The  arguments  of  SETUP  include  control  options  that  describe  the  geometry  as 
well  as  some  computer  operation  parameters.  In  addition,  the  user  must  also  provide  a 
subroutine  GETUM  for  data  input. 
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MAIN  PROGRAM 

PARAMETER  (NDIMX=  180) 

DIMENSION  X(NDIMX,NDIMX) 

COMMON  /GET/KDUM 
COMMON  WORK  (300000),ANG(50) 
COMMON/OUTCOM/LUNOUT,I80132 
COMMON/PARM/IPAR(I2),PAR(3) 

EQUIVALENCE  (NDIMU,IPAR{l)),(ICIR,IPAR(2)),(IGEOM,IPAR(3)). 

1  (NANG,IPAR(4)),(MODANG,IPAR(5)),(KDIMU,  IPAR(6)). 

2  (IMIT,IPAR(7)),(NWORK,  IPAR(8)).(NFLOAT,  IPAR(9)), 

3  {ISTORE,IPAR(IO)).(IPRINT,IPAR(1 1)),(LUNATN,IPAR(I2)), 

4  (PWID,  PAR(1)),(AXISU,  PAR(2)),(RFAN.  PAR(3)) 

EXTERNAL  PRFF.BRFF 

LUNOUT  =  6 

180132  =  1 

NDIMU  =  180 

ICIR  =  1 

IMIT  =  0 

IGEOM  =  2 

KDIMU  =  235 

NANG  =  15 

MODANG  =  0 

NWORK  =  300000 

LUNATN  =  30 

NFLOAT  =  1 

ISTORE  =  0 

IPRINT  =  4 

RFAN  =  4367.3364 

PWID  =  1.3 

KDUM  =  KDIMU 

DIMUl  =  KDIMU 

AXISU  =  (DIMUl +  l.)/2. 

CALL  SETUP(IPAR,PAR,ANG) 

ISTP  =  10 
IRLX  =  1 
lERR  =  0 
IZER  =  0 

CALL  CONGR(X,PRFF,BRFF,ISTP,IRLX,IERR.IZER) 

22  STOP  23 
END 
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X 

PRFF 

BRFF 

LUNOUT 

180132 


NDIMU 

ICIR 

IMIT 

IGEOM 

0 

1 

2 

3 

KDIMU 

NANG 

MODANG 

0 

1 

2 

3 

4 

5 
-I 


NWORK 

LUNATN 

NFLOAT 


PARAMETER  DEFINITION 

THE  RECONSTRUCTION  ARRAY 
THE  PROJECTION  SUBROUTINE 
THE  BACK-PROJECTION  SUBROUTINE 
LOGICAL  UNIT  NUMBER  FOR  OUTPUT 

FLAG  INDICATING  NUMBER  OF  CHARACTERS  IN  A  LINE  OF 
OUTPUT  ON  LUNOUT 

0  =  80  CHARACTERS  (132  CHARACTERS  OTHERWISE) 

THE  LINEAR  DIMENSION  OF  THE  RECONSTRUCTION  ARRAY 
0  TO  RECONSTRUCT  A  CIRCULAR  ARRAY  OTHERWISE 
RECONSTRUCT  A  SQUARE  ARRAY 

0  TO  RECONSTRUCT  EMISSION  DATA  OTHERWISE 

RECONSTRUCT  TRANSMISSION  DATA 

GEOMETRY  FLAG 

PARALLEL  BEAM  GEOMETRY 

FAN  BEAM  GEOMETRY  (CURVED  DETECTOR) 

FAN  BEAM  GEOMETRY  (FLAT  DETECTOR) 

RING  DETECTOR  GEOMETRY 

NUMBER  OF  BINS  IN  THE  PROJECTION  ARRAY  SUPPLIED  BY 
THE  USER 

NUMBER  OF  PROJECTION  ANGLES 

MODE  FOR  PROJECTION  ANGLE  INPUT 

USER  SUPPLIES  PROJECTION  ANGLES  IN  DEGREES 

USER  SUPPLIES  PROJECTION  ANGLES  IN  RADIANS 

PROJECTION  ANGLES  GENERATED  BETWEEN  ZERO  AND  PI 

STARTING  AT  THE  HALF  ANGLE 

PROJECTION  ANGLES  GENERATED  BETWEEN  ZERO  AND  2*PI 
STARTING  AT  THE  HALF  ANGLE 

PROJECTION  ANGLES  GENERATED  BETWEEN  ZERO  AND  PI 
STARTING  AT  ZERO 

PROJECTION  ANGLES  GENERATED  BETWEEN  ZERO  AND  2*PI 
STARTING  AT  ZERO 

WHERE  I  IS  BETWEEN  2  AND  5  DOES  THE  SAME  AS  ABOVE  WITH 
THE  ORDER  OF  ANGLES  REVERSED 

DIMENSION  OF  THE  USER’S  COMMON  BLOCK  IN  BLANK 
COMMON 

LOGICAL  UNIT  NUMBER  FOR  ATTENUATION  FACTOR  STORAGE 
NUMBER  OF  WORDS  FOR  A  FLOATING  POINT  VARIABLE 
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ISTORE  -  0  TO  PERFORM  A  RECONSTRUCTION  OTHERWISE  ONLY  DO  A 

STORAGE  SIZE  TEST 

WORK  -  BLANK  COMMON  WORKING  ARRAY 

IPRINT  -  PRINT  FLAGS  (BIT  0  =  LEAST  SIGNIFICANT  BIT) 

BIT  0  PRINT  REQUIRED  FLOATING  POINT  BLANK  COMMON 
WHENEVER  CHANGED 

BIT  1  PRINT  PROJECTION  DATA  AND  UNCERTAINTIES 
BIT  2  PRINT  SETUP  VALUES  FROM  IPAR  AND  PAR  ARRAYS 
BIT  3  PRINT  FILTER  FUNCTION  FOR  CONVOLUTION  AND  FILTER 
ROUTINES 

BIT  4  PRINT  VALUES  FOR  THE  LAGRANGE  MULTIPLIERS  AND 
THE  GRADIENT  FOR  THE  FUNCTION  OF  LAGRANGE 
MULTIPLIERS  FOR  THE  ENTROPY  RECONSTRUCTION 
BIT  5  PRINT  POINTERS  IN  BLANK  COMMON 
RFAN  -  FOR  FAN  BEAM  GEOMETRY  RFAN  IS  THE  DISTANCE  FROM  THE 
SOURCE  TO  THE  CENTER  OF  ROTATION.  RFAN  IS  MEASURED 
IN  UNITS  OF  PROJECTION  BIN  WIDTHS  AT  THE  CENTER  OF 
ROTATION. 

PWID  -  PIXEL  WIDTH  (IN  UNITS  OF  PROJECTION  BIN  WIDTH) 

AXISU  -  THE  PROJECTED  LOCATION  OF  THE  ROTATION  AXIS  IN  THE 
PROJECTION  ARRAY  (THIS  IS  SUPPLIED  BY  THE  USER  AND  IF 
AXISU  IS  INTEGER,  THEN  ROTATION  AXIS  FALLS  IN  THE 
CENTER  OF  A  PROJECTION  BIN.) 

ISTP  -  NUMBER  OF  ITERATION  STEPS 

IRLX  -  IRLX  IS  NON-ZERO  FOR  ITERATIVE  RELAXATION 

lERR  -  lERR  IS  NON-ZERO  FOR  WEIGHTED  LEAST  SQUARE 

IZER  -  IZER  IS  EQUAL  TO  0  IF  INITIAL  SOLUTION  EQUALS  0 

IPAR  -  INTEGER  ARRAY  WHICH  CONTAINS  THE  OPTIONS  AND 

VARIABLES  DESCRIBED  IN  THE  EQUIVALENCE  STATEMENT 
PAR  -  REAL  ARRAY  WHICH  CONTAINS  THE  PARAMETERS  DESCRIBED 

IN  THE  EQUIVALENCE  STATEMENT 
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